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In this work we extend the etching model [T] to d + 1 dimensions. This permits us to investigate 
its exponents behaviour on higher dimensions, to try to verify the existence of an upper critical 
dimension for the KPZ equations, with our results sugesting that d = 4 is not an upper critical 
dimension for the etching model. 


I. INTRODUCTION 

Among the many equations describing the behaviour 
of moving surfaces, the KPZ equation [5] is of special 
interest, partially motivated by the fact that it describes 
many physical phenomena, such as flame front propaga¬ 
tion laii! and deposition of thin films [5]. Although its 
exponents in one dimension are long known 0 |7j, the 
renormalization group technique used for obtaining its 
values is not usable for higher dimensions. As such, usu- 
aly either numerical simulations or approximate meth¬ 
ods are used to obtain exponents for the KPZ in higher 
dimensions, and even then numerical determination of 
exponents can be tricky [8]. 

Such methods give no final answer to the values of 
these exponents, which combined with the absence of ex¬ 
act solutions leads to the much debated possibility of 
an upper critical dimension (UCD) dc for the dynamic 
exponents [5] . 

One approach consisted of trying to find a general ex¬ 
pression for the exponents depending on the substrate 
dimension d, i.e. a = a{d), /3 = fi{d) and z = z{d). No¬ 
table examples are those for the RSOS model, by Kim 
and Kosterlitz nni, for the Eden model El by Kertesz 
and Wolf, the heuristic approach to the strong-coupling 
regime by Stepanow [T^, a tentative method based on 
quantization of the exponents by Lassig[T3] and a pertu- 
bation expansion of the KPZ equation by Bouchaud and 
Cates [H]. Unfortunately, further numerical results have 
shown these results to be incorrect [IMIH]- 

Analytical methods such as mapping of the directed 
polymer [19], perturbation expansion |20| and mode¬ 
coupling techniques [2T] among others [22U25] observe 
dc = 4. A asymptotic weak noise approach by Fogedby 
[5S| suggests dc < 4. Otherwise, some numerical studies 
found no such limit dUEiiES], as well as the numerical 
and theoretical results by Scharwartz and Perlsman |53| . 

We present in this work a version of the Etching model 
by Mello et al [H ED] extended to work with d -I- 1 spa¬ 
tial dimensions |31| . The one dimensional version of the 
Etching model is known to be compatible with the KPZ 
equation, and as such, is classified in the KPZ universal- 
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ity class. 

Using this version of the model, we determine expo¬ 
nents for surfaces with 1 < d < 6 , reaching the conclu¬ 
sion that if there is a UCD, it must be such that dc > 6. 


II. THE ETCHING MODEL IN d DIMENSIONS 

Proposed by Mello et al [T] in 2001, the etching model 
simulates a one dimensional crystalline solid submerged 
in a solvent liquid. Its scaling exponents are very close to 
those of the KPZ equation, namely a = 0.4961 ± 0.0003 
and /3 = 0.330 ± 0.001, and as such, it is believed do 
belong to the KPZ universality class. This model was 
object of extensive research in recent years [32H37] . 

We extend the model to d -I- 1 dimensions, considering 
the solid a square lattice exposed to the solvent liquid, 
with a removal probability proportional to the exposed 
area. The algorithm can be described as 

1. at discrete instant T one horizontal site i = 1, 2..., L 
is randomly chosen; 

2. /i,(r+l) = /i,(T) + l; 

3. if hi+s{T) < hi{T), do hi+s{T -f 1) = hi{T), where 
(5 = ±1 are the first neighbours. 

In the multidimensional version i and S are vectors and 
6 runs over the 2‘^ first neighbours of the hypercube. We 
consider L to be the substrate length in each direction, 
with the total number of sites is L‘^. The normalized 
time t defines the time unity as L‘^ cellular automata 
iterations, i.e., t = T/L'^. We use periodic boundary 
conditions on the surface to reduce unwanted finite size 
length effects. Albeit the model is not a direct mapping 
of the KPZ equation, it generally mimics its dynamics 
and reproduces its exponents for 1-1-1 as well as the 
general case d -I-1. 

We simulate several substrate lengths L for each di¬ 
mension d, with each experiment being repeated sev¬ 
eral times. This ensemble average is necessary to reduce 
noise, producing higher accuracy in the resulting expo¬ 
nents. 

On figure]^ we show our results for simulating rough¬ 
ness dynamics on various substrate lengths for dimen¬ 
sions from 1 -I- 1 to 6 -I- 1, on log-log scale. On all simu- 
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lated dimensions, the expected Family-Vicsek (FV) scall- 
ing [38] is visible. 

The Family-Vicsek scaling is a relation that can be 
used to model surface rougthness dynamics by consider¬ 
ing it composed of two different regimes: one in which 
it grows in a power-law like function of time , and af¬ 
ter a saturation time tx) h saturates. The values of the 
saturation roughness is related to the substrate length, 
with 

Ws oc L“, ( 1 ) 

These properties are expressed in the FV relation: 

/ ../ / I if t tx 

w{t,L) = Wsf{t/t-^,P) = I y + 

I Ids it r » Ex 


III. DETERMINATION OF EXPONENTS 
VALUE 

Using the FV relation, we obtain our exponents by 
fitting our data to a set of power laws. We fit values of 
Wy, /3l, and Wg by using the two expressions of ® at 
t tx and t ^ tx- Determination of tx is made by 
analysing the intersection between the functions of the 
aforementioned regimes. 

The etching model presents a transient time at the 
beginning of the growth process, and as such data from 
these times where t < 1 are discarded from the fitting. 
This implies in a /3 that is not independent of L. For this 
reason, the parameter obtained from the fitting is called 
pL and a correction is made, in the form: 

/3l = ^ , (3) 

where 7 « 1 . 

This correction considers the real value of /3 to be the 
asymptotic of /3l, eliminating finite size effects. 

The values of /3 l, Wg and tx for each value of d and 
L were obtained from roughness fitting and plotted in 
figure 

IV. DYNAMIC EXPONENTS AND THE UCD 

Using the fitting from the data shown on figure we 
obtain values for each exponent for dimensions ranging 
from 1 -I- 1 to 6 -I- 1. It allows us to observe how these 
expoents behave on higher dimensions. 

On table |T| we show our results. It is simple to observe 
that the expected behaviour for a system with a UCD 
on dc = 4 is nowhere to be found, with all exponents 
continualy changing. This suggests that there is not an 
UCD dc < 6 , in agreement with previous results [niiiTi- 
Uni I2S| obtained through other models. 
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FIG. 2: Parameters of ([^ plotted as functions of L with 
d = 1... 6. (a) Ws{L), (b) tx (L) multiplied by 10“* for better 
visualization, and (c) Pl{L). In all simulations L = 2", with 
n integer. 

TABLE I: Dynamic exponents obtained from the fittings of 
figure A evidence of the precision of these exponents is the 
value of a + z, which should be 2. 

da j3 z a + 2 

1 0.497(5) 0.331(3) 1.50(8) 2.00(1) 

2 0.369(8) 0.244(2) 1.61(5) 1.98(2) 

3 0.280(7) 0.168(1) 1.75(9) 2.03(2) 

4 0.205(3) 0.116(3) 1.81(3) 2.02(1) 

5 0.154(2) 0.079(3) 1.88(6) 2.04(1) 

6 0.117(1) 0.054(1) 1.90(6) 2.01(1) 


It is important to note that although the one dimen¬ 
sional Etching model is on the KPZ universality class, it 
is hard to classify it, or for that matter, any multi di¬ 
mensional model, on the KPZ universality class, as there 
are no known solutions for this case. We can, however, 
compare our results with with others m, showing great 
concordance. 

V. CONCLUSION 

We have made a generalized version of the etching 
model, capable of simulating surfaces on d -|- 1 dimen¬ 
sions. Using this version of the model, we were capable 
of obtaining exponents for systems up to 6-1-1 dimensions. 

Throught the obtained exponents, we have shown that 
the etching model does not show a upper critical dimen¬ 
sion at dc = 4. It is not possible, however, to assert 
that the same thing is true to the KPZ equation, as we 
still do not have a formal mapping of our model to the 
KPZ equation, although comparison with current liter¬ 
ature hints into the same general direction of both the 
etching model belonging to the KPZ universality class 
and the absence of this UCD on the KPZ equation. 

We expect as well that importa results obtained in 
stochastic process [iOHiS] could be used to give a more 
complete solution to this problem. 
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